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Abstract A linear quadratic Dirichlet control problem posed on a possibly non-convex 
polygonal domain is analyzed. Detailed regularity results are provided in classical 
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1 Introduction 


The investigation of optimal control problems with partial differential equations has been 
of increasing interest in the last decades. In this paper we will study the control problem 


mm, 


(P)< 


in J(u) = 7 - [ (Su(x) - yn(x)) 2 dx + ^ [ u 2 (x)da(x) 
4 Jn Dr 

subject to ( Su,u ) £ 7/ 1 / 2 (D) x L 2 (T), 

, u £ U a d = {u £ L 2 (T) : a < u{x) < b for a.a. x £ T}, 
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where Su is the solution y of the state equation 

— Ay = 0 in D, y = u on T, (1.1) 

the domain 11 C M 2 is bounded and polygonal, T is its boundary, a < b and v > 0 are real 
constants and t/q is a function whose precise regularity will be stated when necessary. 
Note that, for u E U a d , the state equation does not possess a variational solution in 
general, so a very weak solution is considered (see Theorem 2|). We will discuss here the 
regularities of the optimal state y, the optimal control u and the corresponding adjoint 
state (p which are limited by singularities due to corners of the domain and due to the 
presence of control constraints. 

The classical Sobolev (Slobodetskii) spaces are denoted by W t,p (D,) and, in the case 
of p = 2 , by 1^(11). As usual, for t > 0 , Wq P (D) or Hq{PI) will denote the closure 
respectively in IT i,p (n) or i^ f (fi) of D(Ll), the space of infinitely differentiable functions 
with compact support in 17, and W~ t,q (D ) with q+ p~ 2 = 1 [resp. i7 _4 (n)] is the dual 
space of Wq V (D) [resp. Hq(D)]. 

The seminal paper on Dirichlet control problems is the work by Casas and Raymond 
[ 6 ]. They investigate the problem even with a semilinear state equation. Assuming a con¬ 
vex polygonal domain with maximal interior angle uq < 7 r, they prove u E kF 1 ~ 1//p,p (T) 
and y E W 1,p (Pl) withp < pn = 2/(2 — min{Ai, 2}) and Ai = n/ui both for the case with 
and without control constraints. Note that pn > 2 due to the convexity of the domain. 
May, Rannacher and Vexler [14] consider the unconstrained Dirichlet control problem 
with linear state equation m also in convex domains and derive u E H 1 1 //p (r) and 
y E -£f 3 / 2 - 1 /P(f7) for p < pq. Deckelnick, Gunter and Hinze [9] focus on approximation 
issues in the case of smooth domains (class C 3 ) in two and three space dimensions. The 
regularity is determined by the box constraints since corner (or edge) singularities do 
not occur. Finally we would like to mention the paper m by Of, Phan, and Steinbach 
where the control is searched in such that the state y E i/ 1 (D) satisfies the 

weak formulation; in this case the regularity issues are less severe. 

Due to the assumptions on the domain, the publications El, IS], and m have in com¬ 
mon that the adjoint problem can be solved in H 2 (Ll) Pi Hq (O). For that reason the very 
weak formulation of the state equation is well defined in these publications with test 
functions from H 2 (Q) I ~1 Hq(D). This is not the case when non-convex domains are con¬ 
sidered. Instead, the very weak fomulation of the state equation should be defined with 
test functions from H^(Q) n Hq(Q), where H^(Q) := {u E I/ 1 (D) : Av E T 2 (0)}, see 
the paper by Casas, Mateos and Raymond [5j (A. 16)] and also [lj for further approaches 
how to understand the solution of the Poisson equation with non-smooth boundary data. 

In the paper at hand we prove basic regularity results for the solution of the state and 
adjoint state equations in Section[2] We extend in Theorem[I]to non-convex domains the 
well-known H l / 2 (D) regularity of the solution of the state equation (11.11) and we prove 
in Theorem [7] that the maximum principle also holds for very weak solutions. Our main 
regularity results for the variables in the optimal control problem are proved in Section 
[3] Among many other results, we prove that, in the presence of control constraints and 
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under minimal regularity assumptions on the data (jjq £ L s *(ft), s* > 2), the optimal 
control is continuous despite the possible non-convexity of the domain (cf. Theorem 
ED. The main idea is very simple: we compute explicitly the normal derivative of the 
singular part of the adjoint state and exploit the projection relation established by the 
first order optimality conditions. We also investigate the case s* = 2. This case has not 
been treated by the cited references for convex domains. In Section 0] we prove that, for 
regular data, the regularity of the optimal solution is indeed slightly better. We give 
conditions for the control to be in Pf 3 / 2_e (ft) for all e > 0 (Corollary [TUD, which is the 
best regularity we can expect under pointwise control constraints. These results will 
be helpful to derive error estimates for finite element approximations of problem (P). 
The numerical analysis will be carried out in a forthcoming paper. The short Section 
[5] is devoted to the unconstrained case. Whereas in many other control problems the 
regularity of the unconstrained solution is better than that of the constrained solution, we 
encounter here the phenomenon that the constraint inhibits poles of the unconstrained 
solution. For that reason the regularity of the constrained control is determined, up to 
some exceptional cases (cf. Remark 1121) . by the largest convex angle but the regularity 
of the unconstrained control is determined by the overall largest angle (cf. Corollary 

n. 

Finally we would like to remark that sometimes the state equation is considered in 
the form 

- Ay = / in ft, y = u + g on T, (1.2) 

but we can take the unique function yo which solves — Ayo = / in ft, yo = 9 on T, and 
replace y := y + yo, yn ■= yn ~ yo and recover problem (P) with (11.11) for data sufficiently 
smooth. 


2 Notation and basic results for elliptic equations 

Let us denote by M the number of sides of T and {xj ( its vertexes, ordered coun¬ 
terclockwise. For convenience denote also xq = xm and xm+i = x\. We will denote 
by Fj the side of F connecting Xj and Xj + \, and by uij £ (0, 27 r) the angle interior to ft 
at Xj, i.e., the angle defined by Fj and r.,_i, measured counterclockwise. Notice that 
To = Fm- We will use ( rj,9j ) as local polar coordinates at Xj, with rj = \x — Xj\ and 
9j the angle defined by Tj and the segment [xj , x]. In order to describe the regularity of 
the functions near the corners, we will introduce for every j = 1 ,,M the infinite cone 


Kj = {x £ M 2 : 0 < rj, 0 <6j < ojj} 
and a positive number Rj such that the sets 

Nj = {x £ M 2 : 0 < rj < 2Rj, 0 < 9j < ujj}, 
satisfy Nj C ft for all j and IVj fl Nj = 0 if i ^ j. 
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For every j = 1 ,,M we will also consider 


: ► [ 0 ; 1 ] 


an infinitely differentiable cut-off function which is equal to 1 in the set {x G M 2 : rj < 
Rj} and equal to 0 in the set {x G M 2 : rj > 2 Rj}. 

We will denote by Zf the solution of the homogeneous Dirichlet problem with dis¬ 
tributed data 

— Az = f in 12, z = 0on T. (2-1) 

The regularity of / and Zf , as well as in what sense the equation must be understood, 
will be stated when necessary. 

For every j = 1 ,,M we will call A j the leading singular exponent associated with 
the operator corresponding to the corner Xj. For the Laplace operator it is well known 
that Xj = n/ujj. For convenience we will suppose that Ai = min{Aj : j = 1,..., M}. In 
general, the maximum regularity in non-weighted Sobolev spaces of the solution of the 
Poisson problem for regular data will be given by the Sobolev exponents 


Pn 


2 

2 — min{Ai, 2} 


and tn 


1 + Ai. 


This is, for / G C°°(D), the solution zj of equation (12.11) will satisfy that z € W 2,p (0) D 
if*(fI), for all p < pn (cf. Grisvard [Till Th. 4.4.3.7]) and all t < tn (cf. Dauge 0 
§23.C]; see also Grisvard [11] Th. 2.4.3 and §2.7]). For less regular data / G W~ 1,q (D), 
the maximum regularity of the solution of the homogeneous Dirichlet problem is given 
by the exponent 

2 

PD = -j- ■ r, ' -i • 

1 — mmjl, Ai) 

This means that Zf G W l,q {0) if q < pd (cf. Dauge [8] Theorem l.l(i)]; see also Jerison 
and Kenig [121 Thms. 0.5, 1.1, 1.3]). Further results on the local regularity in each 
corner can be stated using weighted Sobolev spaces. 

We will state now two lemmas collecting several regularity results that will be used 
later to state the regularity of the solution of the control problem. The following lemma 
is a consequence of well known regularity results collected in the book by Grisvard ma¬ 
lt tells us accurately how the singularities arising from the corners behave for problems 
with regular data. We introduce the following sets in order to describe the singular 
behavior of the solution of the Poisson problem at the corners. 

For 1 < p < +oo such that 


and m G Z, define 

J™ 


2(P-1) 

pXj 


£ZV?G M} 


l j G {1,... , M} such that 0 < m\j < 2 



( 2 . 2 ) 


(2.3) 
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The condition in (12.211 is necessary to deduce the TT 2,P (^-regularity of the regular part 
of the solution in the lemma below. The meaning of the sets J1™ is the following: 

£jr™ Xi <£ W 2 ’ P (D) for all j G I™. 

Notice that, for all p < +oo, J™ = 0 for m > 4 since A j > 1/2 for every possible angle 
and 2 — 2/p < 2 for all p < +oo. We also remark that C J 2 C 

Lemma 1. Consider 1 < p < +oo satisfying (12.211 and f G L p (Ll). Then there exist 
unique real numbers (cj >m )j e j™. and a unique solution Zf G H 1 (f2) of problem (12.111 such 
that 

3 

Zf = Z reg + Cj tm £,jr™ 3 Sm(m,\j8j) 

rn=ljeJ™ 

where z reg G W 2,p {Ll) and the fj are the cut-off functions introduced above. 

Proof. The result is a direct consequence of uni Theorem 4.4.3.7]. This result can be 
applied since Zf G i7 1 (fl) thanks to [TOj Lemma 4.4.3.1]. □ 

Corollary 2. If f G L 2 {fl), then Zf G i7 4 (n) for all t < 1 + Ai and d n Zf G L 2 (T). 

Proof. To prove this fact, we apply Lemma|TJ the regular part is in H 2 (Q); on the other 
hand, since A j > X± > 1/2 and m > 1, then £,jr™ Xj G H t (Ll) for all t < 1 + Ai. From this 
we obtain that Zf G Since we can choose t > 3/2 we have G L 2 (T). □ 

The next result states the regularity of the solution of problems with boundary data 
in W l l ^ v * ,p * (T) with p* > 2. Let us recall (cf. [10, Theorem 1.5.2.3]) that the trace 
of any function z G VF 1 ,p *(fi) is in ,p * (T), the trace mapping is onto, and, for 

p* > 2 , this space can be characterized as 

{ M 

g G lF 1 _ 1 //p *’ p *(Tj) : g is continuous at every corner Xj 
i =l 

For p* = 2 the continuity requirement in the corners can be weakened to an integral 
condition, see [TO, Theorem 1.5.2.3(c)]. 

Lemma 3. Let g G ,p * (T) for some p* > 2. Then there exists a unique solution 

z G W 1,p (Q), for all p <p*, p < pd, of the equation 

—A 2 = 0 in Q, z = g onT. 

Proof. Due to the trace theorem, there exists a function G G VF 1 ,P *(D) such that its 
trace is 7 G = g. Moreover, we have that AG G VF _ 1 ,P *(D). If we define (, = z — G, then 
it satisfies the boundary value problem 

—A£ = —AG in D, ( = 0 on T. 
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Let p < p*, p < pd- Then £ £ WQ 1,p (f2) (see above or cf. Dauge [U Theorem l.l(i)]; see 
also Jerison and Kenig [T2; Thms. 0.5, 1.1, 1.3]) and hence so does z £ W l ’ p (0). □ 

Since the space of controls is L 2 (T), the state equation must be understood in the 
transposition sense. Following P ) IS], for u £ L 2 (T), we will say that y £ L 2 {PL) is a 
solution of (11.11) if for every / £ F 2 (Pl) 

[ yfdx = - [ ud n Zfdcr(x), (2.4) 

Jn Jr 

where Zf is defined in Lemma [I] The definition makes sense thanks to Corollary [5J 

Existence, uniqueness and regularity of the solution in y £ iL 1//2 (fi) if D, is convex 
domains is a classical result and can be proved via transposition and interpolation. Let 
us briefly recall how this result is obtained for a convex domain. Consider the solution 
operator S of (11.11) with Su = y. Due to the Lemma [3] it is 

S £ ^(if 1 / 2 ^),# 1 ^)). (2.5) 

Using the classical transposition method (cf. [13]) we also have that 

S£C(H~ l / 2 ( T),L 2 (D)). (2.6) 

The final result is obtained by interpolation using that 

L 2 (T) = [iL 1/2 (r),iL- 1/2 (r) ] 1/2 and H l /\D) = [H\0), L 2 (D)] 1/2 , 

see e.g. m Chap. 1, Eq. (2.41)]) for the first result and notice that interpolation 
results are valid for spaces posed on Lipschitz domains (cf. [3] Theorem 12.2.7] or |12|). 
We cannot use this scheme straightforward because (|2.6D uses explicitly that for every 
/ £ L 2 (D), Zf £ H 2 (D), and this is not true for non-convex polygonal domains. 

For problems posed on non-convex polygonal domains, Berggren [2j Theorem 4.2] 
proves existence, uniqueness and regularity of the solution y £ 7L f (D) for every 0 < t < 
e < 1/2 (e depends on the domain). We can also achieve y £ L7 1//2 (D) in non-convex 
domains. The proof of our following result uses interpolation spaces; a different proof 
by using integral operators is given in JT]. 

Theorem 4. For every u £ L 2 (T) there exists a unique solution y £ H l ^ 2 {D) of (11.11) 
and 

llylliU/2(n) < C1MlL 2 (r)- 

Proof. Notice that for 0 < e < 1/2, we also have that, for 6 = 1/(1 + 2e), 

L 2 (r) = [H 1/2 (T),H~ £ (T)} g and H ^(D) = H 1/2 ~ £ (Q)] e 

and therefore the result will be true if we can prove that S £ £(H~ £ (F), H l / 2 ~ £ (Ll)) for 
some e > 0. 
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Fix 0 < e < min{Ai — 
every / € L 2 {Pl) 


1/2,1/2}. For any u £ H £ (r), we will say that y 



-(u,d n Zf) H -e(T) H e(Yy 


Su if for 
(2.7) 


Notice that since e < Ai—1/2 we have that Zf £ H ?> / 2+e (D) and that d n zj £ FI /= i = 

H £ (T) because e < 1/2 (cf. [10, Theorem 1.5.2.3(a)]). Moreover, due to the trace theo¬ 
rem and elliptic regularity, we have that 


\\9 n Zf\\HRr) < C\\f\\ H e-i/ 2 {p,y 


( 2 . 8 ) 


Notice also that if u £ H l D(T) then the unique variational solution y £ N 1 (Q) of 


— A y = 0 in 0, y = u on T 


is a solution of (12.711 and if u £ L 2 (T), then (12.71) is the same as (12.41) . 

Let us prove uniqueness of the solution of (12.71) in L 2 (Q) first. If u = 0 and y = Su £ 
L 2 (Q), then, taking / = y as test function in (12.71) we get f n y 2 dx = 0, and therefore 
y = 0. Since the problem is linear, the solution is unique. 

We next prove existence of a solution y £ H l / 2 ~ £ (Ll) of (12.71) . We know (cf. [TO , 
Theorem 1.4.2.4]) that H l / 2 ~ £ {D) = iT 0 1/2 “ £ (D) and hence (i7 1 / 2 - £ (0))' = H £ ~ l / 2 {D). 
Denote 

F ={f £ L 2 (D) such that ||/|| He -i/2 (n) = !}• 

For any u £ H l / 2 (T ), y = Su, we use that L 2 (D) is dense in H £ ~ l / 2 {Pl) and (12.81) to 
obtain 


IMI H i/ 2 - e(Q) = sup(/,y) HS -i/2 (n)) Hi/a-e ( n ) = sup / yfdx 

j j «/ Q 

= sup -(u,d n Zf) H -e(ry H e(r) < sup ||<9n^|i^(r)IMItf-*(r) 

< C sup ||/|| ff£ -l/ 2 (Q) = C'||u||j i --e( r ). 

f 


(2.9) 


Take a sequence u k £ H 1 I 2 ( T), u k —>■ u £ H £ (T), and let y k = Su k . We have just 
proved that 

\\Uk — ym||/fi/ 2 - e (f2) — ll^fc — u m\\H- s (T) 

and therefore y k converges in H l l 2 ^ £ (PL) to some y £ H 1 / 2 £ (D) that is (the unique) 
solution of the equation: 



/ hill ykfdx = lirri / y k fdx 
In k k Jn 


lim -(u k ,d n z f ) H -e {r y H e {r) = -{limu k ,d n z f ) H -e^ ))H e {r) 
— {u, d n Zf) H -e(Y),H e (r)- 


Finally, (12.91) implies that S £ C(H £ (T),H 1 I 2 £ (D)) and the proof is complete. 


□ 
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The next result is rather technical. It will be used to describe precisely the structure 
of the optimal state in the proof of Theorem [T3] below. In that result we will be able to 
write the control as the sum of a regular part and a singular part. We show in Lemma 
2.5 how to solve the state equation for singular boundary data. Besides the usual regular 
and singular parts that we described in Lemma [1] a new singular part of the solution 
arises from the boundary data. 

Define the jump functions at the corners 


f 1 on {x G dKj : 6j = 0} 

( —1 on {x G dKj : 0j = Uj}. 


( 2 . 10 ) 


Notice that Xj = 1 on and Xj = — 1 on Tj_i. 

Lemma 5. Consider any pair of subsets H 1 , H 2 C {1 ,,M} and real numbers —1/2 < 
r)j } i and <i 3 .\ for all j G H 1 and 0 < 77^2 and 0^2 for all j G H 2 such that Vj,n/^j ^ ^ for 
any j G BP 1 , n = 1,2. Define 


11 X °i.i ^ r ?’ 1 + X \j a .i-Ai r j J ' 2 on v - 

jeH 1 jew 2 


take p such that 


1 < p < inf 


1 - Vp 


■ j€ 


n = 1,2 and rjj }n < 1 > , 


where we consider inf 0 = +00 and define J™ as in m- Then there exist unique real 
numbers rn = 1,2,3 and a unique solution y G H l / 2 (D) of equation (11.11) 

such that 


y — Vreg + X X 

n=ljgH" 


.Nrn 
*3,^3' j 


a j,n£j 


*3,n\ v 3 


(% 


+ EE 

m=l jes™ 


-3,m^3 


r™ Xj sm(m\j0j) 


( 2 . 11 ) 


where 

Sj, n {0j) = ^ ~ Sin ^ j9j ^ + C ° S ( 2 - 12 ) 

and y reg G W 2 ' p (Pl). If, further, rij^ > 0 for all j Gi 1 , then y G iL 1 (D). 

Proof. Since 77^1 > —1/2, we have that u G L 2 {P) and y G H l / 2 (D) thanks to Theorem 
HI If also r)jx > 0, then u G H 1 /, 2 (T) and hence Lemma [3] gives us that y G H l (D). 

Notice next that rjj tn > —1/2 implies that 2/(1 — r]j >n ) > 4/3 and p is well defined. 

A direct computation shows that for n = 1,2, yj :7l = r^ J,n Sj %n {9) are respectively the 
solutions of the problems 

-At/j,! = 0 in Kj, y jA = r^' 1 on dKj, 
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-Ay /,2 = 0 in Kj, y ji2 = Xjf] 3 ' 2 on 0K r 
Since Ay/,i = 0, we have that 

2 n 

A EE £jVj,n = EE T 2Vy/ )7 jV£/). 

n= 1 / G J n J = 1 j Si 71 


Since |Vy/, n | < Cr^ 3 ’" \ the condition imposed on p implies that 

2 

/ = a EE 

n=l /el” 

and we can write y = y/ + Yll=i Z)/eH" a j,nf.jyj,n , where 

-Ay/ = / in Q, y f = 0 on T. 

Applying Lemma Q] we obtain that y/ = y reg + Y^m=\ J2je jr™ sin(mA/0/) and 

the proof is complete. □ 

Although the maximum principle is a well known result for weak solutions of equation 
CEO (see the celebrated paper by Stampacchia ESI). we have not been able to find a 
reference of its validity for solutions defined in the transposition sense (|2.4|) . For the 
sake of completeness, we include such a result. First, we prove the following technical 
lemma. 

Lemma 6. Consider f £ L 2 (fl), / > 0 and let Zf £ iLg(fl) be the solution of equation 
m- Then d n Zf < 0 a.e. onT. 

Proof. Take u £ C°°(r), u > 0. Thanks to Lemma [3J the solution of equation (11.111 satis¬ 
fies y £ iL 1 (n) and the maximum principle for weak solutions, as proved by Stampacchia 
ESI, holds. Therefore y > 0. Integration by parts shows then that 



and the result follows by the usual density argument. □ 

Theorem 7. Consider u £ L°°(T). Then the solution y of equation (11.11) belongs to 
L°°(Ll) and 

l|y|k°°(Q) < IMIz,°°(r)- 

Proof. Define K = ||u||^oo(pp We will prove that y < K a.e. on D, the proof for —y < K 
being analogous. 
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We already know that y E ^ L 2 (Q) (cf. Theorem 0]). Define Vk = U — K 

and y] t = max(y^, 0) E L 2 (fl). For every / E L 2 (D) we deduce from (12.11) that 


VKfdx 



K)d n Zfda(x ) 


We take / = > 0. For this choice of /, we know from Lemma [6] that <9 n £/ < 0 a.e. 

on r and therefore 

0 < [ (yx) 2 dx = [ yxVxdx = - I (u - K)d n z f dcr(x ) < 0. 

Jn Jq Jr 


So y\ = 0 and y < K. 


□ 


3 Main regularity results for the control problem 

The following result is standard and the proof can be found in [6]. Though in that 
reference only convex domains are taken into account, this result is independent of the 
convexity of the domain, once we have proved Corollary [2] and Theorem 01 Here and in 
the rest of the paper, Proj[ a b ](c) = min{6, max{a, c}} for any real numbers a,b,c. 

Lemma 8. Suppose yn E L 2 (Pl). Then problem (P) has a unique solution u G L 2 (T) 
with related state y G H l i 2 (D) and adjoint state (p G Hq( 0). The following optimality 
system is satisfied: 


u(x) 

= Pr °j [a,b] ( 

—d n ip(x) S j for a.e. x G T, 

(3.1) 

-Ay 

= 0 in LI, y 

= u on T, 

(3.2) 

—A ip 

= y-yn in 

H, (p = 0 on r. 

(3.3) 


As in the proof of Theorem 01 the solution of the state equation (13.21) must be under¬ 
stood in the transposition sense, whereas the adjoint state equation (13.31) has a variational 
solution. 

Next we state a regularity result for the adjoint state in the framework of classical 
Sobolev-Slobodetskii spaces. In the rest of this section we will suppose that yn G L s * (H) 
where, 2 < s* < +oo satisfies 

2(5 *~ 1} gZVjG{l,...,M}. (3.4) 

AjS 

Theorem 9. There exist a unique function <p reg G kF 2,s *(H) and unique real numbers 
where is defined in (12.31) . such that 

3 

= <Preg + ^ ° Sm(m\j9j). (3.5) 

m=i 
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Proof. Since the problem is control constrained, u £ L°°(T), and by the maximum 
principle proved in Theorem [7l y £ Therefore, y — y^ £ L s (li), and we can use 

Lemma [TJ since the related adjoint state (p is the solution of (13.31) we have that there 
exist unique (p reg £ W 2,s (fi) and such that relation (|3.5D holds. □ 

For any s>2we define the set 

= {j £ : A j > 1} 

and for s > 2 and m = 2,3 

Hr = {j € Jr : c j}l = 0}, 

where the coefficients Cj tTn are the coefficients obtained in Theorem [9j Notice that the 
indexes in H), correspond to convex corners. The indexes in Hr correspond to those 
non-convex corners where the main part of the singularity of the adjoint state vanishes, 
and hence the behavior of the solution at those non-convex corners can be somehow 
compared to the behavior of the solution at the convex corners. Notice also that 

U^ =1 ir = U^ =1 Hr (J{j £ I* : Aj < 1 and c jA ^ 0} 

Consider also p > 2 such that for m = 1, 2,3 

P<s*, p< 2 _ 2 mA if j € Hr- (3.6) 

This condition on p appears in a natural way in the proof of Theorem [TT] see (13.101) . If 
s * > 2, then we can choose p > 2. With this choice we have that 

Lemma 10. Let p satisfy (13.61) and for m = 1,2,3 consider j £ H™- Then 

€ W /1 ' 1/p,P ( r )- (3.7) 

Proof. Take j £ U^ =1 Hr and consider Nj the bounded cone of radius 2 Rj defined in 
Section [2j We first prove that Uj = r^ 3 1 £ W ],p (Nj). 

Since j £ U^ =1 Hr, then m\j > 1, so m\j — 1 > 0 and r™ Aj 1 6 C(Nj) C L p (Nj). 

On the other hand iVu^j = (m\j — 1 )r™^ j 2 , and making the usual change of variables 
to polar coordinates, we have 

n n /* 2 J? ' 

J J \Vuj\ p dx = (mXj — 1 ) p ujj J rjr 2 ' )p dr, 

the last integral being convergent if and only if (mXj —2)p+l > — 1. Taking into account 
that j £ Jir implies mXj < 2, the previous condition is fulfilled if and only if p < 9 , 
which is the assumption. 

Relation m now follows from the smoothness of the cut-off function, the trace 
theorem and the continuity of fjr™ 3 *. □ 
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Theorem 11. Let p satisfy (|3.6D . Then, the optimal control u belongs to VF 1—1 / p,p (r), 
the optimal state y belongs to W l,q (D) for all q < p, q < pr>. In particular, if s* > 2, 
both are continuous functions. 


Proof. We will exploit the projection relation (13.11) and the expression for the adjoint 
state obtained in (|3.5I) . 

Notice first that (p reg £ W 2,s * (fl) and <p reg = 0 on T, so d n p reg £ W l ~ l / s *' s * (T) (cf. 
[5j Lemma A.2] for the case s* = 2 or [4] for the case s* > 2). Moreover, if s* > 2, then 
dnTreg{%j ) = 0 on every corner the normal derivative of the regular part is a continuous 
function on T. 

We are going to compute now the normal derivative of the singular part. For any 
m £ {1,2, 3} and j £ J™, we have dn^r™* 3 sm(mXjQj) £ C°°(r \ {xj}), so we have that 
for every compact set /C C T \ {xj : j £ 

d n p€ W l - l ' s *' s \lQ. 

Near the corners, for ry < Rj we have on r ; (where 6j = 0) that 


d n r- lXj =- dgr™ X] sm(mXj9j) 


1 r, m\j . 


= —mXjr " 1X3 1 cos(mAjO) = —mXjr 


mX-i — 1 


(3.8) 


and on (where 9 ; = u )j) that 


d n r™ X:i sm(mXjdj)£j = —dgr" lXj sm(mXj9j) 


1 r, mA, 


= mX,r 


mXn—1 


'p j 


COS 


= (— l) m m.Xjr^ 1 1 . (3.9) 


Next we will distinguish two cases. 

Case 1: if j £ U^ =1 IHI™ then mXj > 1 and hence the limit of both expressions is zero 


x_. ix J ^-''777 = ^xiiL^* ixxv^xx ^ x. cxxxvx xxv^xxv^ ixi^ xxxxxx u ux uw i xx tApxraoxuxxo 10 ZjV^X w 

0. Noticing (13.71) . we have that the choice of the exponent p made in (13.61) gives 


us 


dn Preg + ^ ^ (-j.wfjl'j^ Sm.(mXj9j) J £ W 1 1/p ’ P (T). (3.10) 

So far, we can deduce that, for every compact set K. C T \ {xj : Xj < 1, c ]: \ / 0} 

d n p £ (3.11) 

Case 2: Now j £ Jj», Xj < 1, and Cj t i / 0. We have (—l) m = —1 and mXj — 1 < 0, 
therefore using expressions (13.81) and (13.91) we have 


lim d n Cj j iry’ sin(Xj9j)fj = — sign^ijoo. 
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If it happens that also j € J 2 * U JJ^*, we have that for m = 2 and m = 3, mXj — 1 > 0 
and again the limit of both (13.811 and (13.911 is zero. So we have that 

3 

lim d n ^2 3 sm(m\j6j)£j = - sign(cj,i)oo. 

3 m= 1 


If s* > 2, also on this corner 

lim dn<Preg(x) = 0, 

X^Xj 

and, trivially 

lim d n p(x) = lim I d n (p reg (x) + d n V] sinfmAj^-)^ 

\ m=l 

= — sign(cj j i)oo. (3.12) 

If s* = 2, as we said at the beginning of the proof, d n (p reg € H 1 ^ 2 ( r), so it needs not be 
even a bounded function. Nevertheless, since the singular part behaves like a negative 
power of rj , this term dominates and we also have that (13.1211 holds. 

As a consequence, there exists p g > 0 such that for x € T with |x — Xj\ < p g either 
Proj[ Qi b] d n (p = a or Proj[ a b ] d n (p = b depending on the sign of Cjp. So the control is flat 
near non-convex corners. This, together with the projection formula (|3.1I) and (13.111) 
implies that the optimal control belongs to lT 1_1//p ’ p (r). Finally, the regularity of the 
optimal state y follows from Lemma El □ 


Remark 12. We would like to remark that the case of having dji = 0 can be seen as 
a “rare” case in practice (although this can happen; see Example 1141 belowL So the 
“normal” case is that = 0 for m = 2,3. In this case, in the choice of p made in (13.61) 
the indexes m = 2,3 are excluded, and hence p will only depend maxjwj : ujj < 7r}, so 
we get the same regularity for the control as that obtained in [6] for convex domains. 


To describe more accurately the regularity of the state and the control, we have to 
introduce some further notation. Consider the coefficients c Jim obtained in Theorem [9] 
and define the coefficients 


^j,m — 


mXj Cj^m 

V 

0 


if 0 G [a, b} 
if 0 ^ [a, b} 


(3.13) 


and the functions 

Sj,m(0j) = ^— . ,, s i n (( m Aj - 1 )0j) + cos ((mXj - 1 )9j) . (3.14) 

sm{{rriAj — ijujj) 

The sets will be used now to describe the singular part of the state and the control: 
if j <G H™, then ijr™* 3 1 0 1F 1,S (0) and ^jr^ 3 1 0 VF 1-1 / s ’ s (r). Regarding (13.161) . we 
also mention that if Xj < 1 — 2/s, then 0 kF 1,s (fl). 
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Theorem 13. Assume further that ab 7 ^ 0. Then there exist a unique u reg £ (r), 

a unique y reg £ PF 1,s (fl), for all s < s*, s < po, and, unique real numbers (cj)\.<x-2/s 
such that 

u(x) = u reg + ^2 a L m ^jr mXj ~ l + X/ 1 (3.15) 

m=l,3jeH™ jm 2 st 


3 

y = dreg + ^ ^ a1 '«.).„■ (/>,) + Os/'/' sin ( A A) (3.16) 

m=ljeH™ Aj<l—2/s 

where the \j are the jump functions at the corners defined in (12.101) . the Sj )m {6) are 
defined in ()3. 141) and the fj are the cut-off functions. 


Proof. From the considerations in the proof of Theorem [FT] we have that far from the 
corners with index j £ U^ =1 EI™, the optimal control is the projection of a function 
that is either regular enough or tends to a signed 00 at one point, so it is clear that 

U reg G 

Next we will check what happens in the neighborhoods of the corners with index 

If 0 0 [a, b] then the control would also be flat in neighborhoods of the corners with 
index j € U^ t=1 EI™ again because the normal derivative of the adjoint state is continuous 
near the corner and 0 at the corner. Then (13.151) holds with aj )Tn = 0. 

If a < 0 < b, then the optimal control will coincide with d n (p/v. Using formulas (13.8|) 
and (13.91) and taking into account the definition of the jump functions on the corners Xj 
cnnn, we have that there exists pj > 0 such that for all x £ T such that \x — Xj\ < p g 


u(x) = -d n p(x) = < 


-d n <Preg(x) - C -Alhl r X i 1 
V V J 


-d n p reg (x) -Xg Cj,22Xj rf ] 1 

V V 


3A,-1 


r • 




if j € 

if j £112* u4, 


(3.17) 


and (|3.15|) holds for a^ m = —mcj ym \j/v. 

Let us finally check (13.161) . We will write y = y\ + y 2 , where 


-Ayi = 0 in D, y 1 = u reg on T 


and 


-A y 2 = 0 in Q, y 2 = ^ ^ a j,rrfjr mX] 1 + ^ XYEpffr 2 ^ 1 on T. 

m=l,3/eH™ /eHj* 
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Using Lemma[3]we have that y\ £ W 1 ,s (fi) for s < s*, s < pr>- From Lemma[5]we have 
that there exist a unique 2 / 2 ,reg £ TU 2 ’ P (U), p defined in (13.61) . and unique real numbers 
c,j, m such that 

3 

2/2 = 2/2 ,reg T ^ ^ ^ ^ ) 

m =1 

3 

+ y, r j-t»\i r ';‘ x ' sin ( m0 j) 

m=l jel™ 

Since p >2 and s < s* < oo, in dimension 2 we have thanks to usual Sobolev’s imbedding 
IU 2 ’ P ( 0 ) W 1 ’®^) that 2 / 2 , reg € W 1 ’*^). 

As we already mentioned, among the terms of the second addend, those which corre¬ 
spond to are not in 1U 1,S (U). Notice that if s* < po, and hence s < po , then there 
could be some terms in fU 1 ,, 5 (U), which we gather in a function y a ,reg € kF 1 ,s (U). 

For the last addend, we have that Zjr™* 3 sm(rnOj) ^ PF 1 ,s (0) iff m\j < 1 — 2/s. 
Since s > 2 , this excludes the case m > 1. We gather all the other terms in a function 

2/6,reg € W 1 ’’^). 

So finally we have that the (|3.16l) holds for y reg = 2/1 + 2/2, reg + 2/a, reg + 2/6, reg € kF 1 ,s (U) 
and Cj = Cj.i. □ 


Let us present now the example announced in Remark [T3J For the example, we want 
to remark that our results are also applicable in curvilinear polygons without many 
changes (see m Th. 5.2.7]). The only thing to take into account is that if the angle ioj 
between two curved arcs is grater than 7 r, then we must impose also s* < uij/(ujj — i r) 
(this is not the case in the following example). 

Example 14. Let ui\ = 3tt/2 and consider the curvilinear polygon D = {x £ M 2 : 
0<r<l, 0 < 0 < wi}, where (r, 0) are the usual polar coordinates. We have that 
W 2 = W 3 = 7 r/ 2 , and hence Ai = 2/3, A 2 = A 3 = 2. Suppose yn £ L°°(Q), so we may 
choose any s* < + 00 . We have jJ, = JJ 2 * = {1} and jf* = 0. We also have p < 3. The 
adjoint state can be written as 


V = Preg + Cl,l6 7 ’ 2/3 Sm ( ) + XlCl^W 


A/3 


Sill 


where p reg £ W 2)S (U) for all s* < + 00 . Take —a = b = 1 for instance. If cpi 7 ^ 0 , then 
= 0 for m = 1,2,3 and hence u = u reg £ W l ~ l / s * ,s * (T). 


Define now 


yn(x) 


1 if 9 < ( jJ \/2 

— 1 if 9 > u>\/2 


(3.18) 


such that the problem is skew-symmetric with respect to the line with 6 = uq/2. The 
skew-symmetry of the data suggests that the solution is skew-symmetric, i. e. that 
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the symmetric contribution with r 2 / 3 sin (|$) vanishes, cgi = 0 (result which we have 
confirmed numerically), and hence 

U = U reg + Xiai,2£l7- 1/3 

so u G W l ~ l /P’P{T) for all p < 3. 

In Theorem [13] we have excluded the cases a = 0 or b = 0. These cases can be treated 
with the same techniques. Nevertheless, many cases may appear depending on which 
of the bounds is zero and the sign of the coefficients of the singular part Cj_ m . As an 
example, we will show how to treat some of these cases. We will discuss first what we 
think is the “generic” case, and then a seemingly more “rare” case. Without loss of 
generality suppose a = 0 , b > 0 . 

Case 1 Take j G H]* and suppose dj } \ ^ 0. Then in the expression for the normal 
derivative of the adjoint state 113.171) . the term £yry' 1 dominates the term d n ip reg , since 
^ ^i-i/s*,s*(r) and d n {p reg G W 1_1 / s * ,s *(r). Therefore, if c 3 ,\ < 0, we would 
have that d n tp(x) G [0,6] in a neighborhood of Xj, and hence u{x) can be computed as 
in ([3.171) . On the other hand, if Cj.\ > 0, then d n (p{x) < 0 in a neighborhood of Xj, so 
u(x) = 0 in that neighborhood. This would be the case of taking Example QT] with the 
following data: = 37 t/4 , a = 0, 6 = 1 and either = 1 or ijq = — 1, which would 

give sign(cyi) = - sign (y n ). 

Case 2 If j G H 2 , and c J: 2 < 0, then d n cp(x) G [0, 6 ] in a neighborhood of Xj on the side 
Tj, but d n (p{x ) < 0 in a neighborhood of Xj on the side Tj_i, so on Tj, u(x) would have 
the same expression as in (13.171) . but, on Tj_i, u{x) would be flat near the corner Xj. 
This would be the case of taking Example [T4l for a = 0, 6 = 1 and y q defined in ([3.181) . 


4 More regular data 

Taking advantage of the regularity of the optimal state, we can obtain several results 
for more regular data. We will write some results that we think will be useful for the 
numerical analysis of problem (P). To be specific, to obtain error estimates for problems 
with regular data, we will need that for y n G H l {D), u G H 3 / 2 ~ £ (T) (cf. Corollary fTHI) 
and VE 3 ,p (fI) regularity of the regular part of the adjoint state if yn G W 1)P {0), p > 2 
(see Corollary [TT]) . Some other results can be obtained used the techniques exposed in 
Section [3] 

Suppose now that yn G H 1 * (Ll), with 0 < t* < 1 such that 

(1 + t*)/Xj ^ Z Vj G {1, ... ,M}. 

For t > — 1 and mGZ define 

J™' = { j G {1 ,,M} such that 0 < m\j < 1 + 1} 
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Notice again that due to our choice of t*, we only will deal with the cases m = 1,2,3. 

Corollary 15. There exist a unique function (p reg £ H 2+t (0) and unique real numbers 
c ]m such that 

3 

T = Treg + ^ ° SW.(m\j9j). 

m=1 


Proof. Since t* > 0, there exists s* > 2 satisfying (13.41) such that yo £ L s * (Q), and hence 
we can apply Theorem 1111 and we have that also y £ W 1,q (D) ^ H 1 (fl) ^ H 1 (0), 
where q > 2 is defined in Theorem El The result follows directly from the adjoint state 
equation (|3.3p thanks to the regularity results in [7J §23.C]; see also JTT, Th. 2.4.3 and 
§2.7]. Notice that due to the conditions imposed on t*, logarithmic terms do not appear 
in the development of the singular part. □ 

To describe the regularity of the optimal control and state, we first introduce the sets 
in an analogous way as we did for the sets H" 1 , the indexes being taken now in the 
sets defined above instead of the sets J” 1 and considering the coefficients c ]t \ obtained 
in Corollary El We next define the exponents t > 0, related to the regularity of the 
control, and t > 0, related to the regularity of the state such that 

t < t*, t < 1, t < m\j — 1 if j £ H" 1 , m = 1,2, 3 
t < t, t < Ai. 

The meaning of these bounds is the following. The regularity of the optimal control will 
be limited by the regularity of the data, the impossibility of having a control globally in 
77 3 /2 (r) ^ ue £j ie corners anc [ bound constraints, and the singular behavior of the 
control at the convex corners or the “special” nonconvex corners that may lay in Wf 1 
for m = 2,3. The regularity of the optimal state will be limited by the regularity of the 
control and the singular behavior at the nonconvex corners of the solution. 

Corollary 16. Suppose that 

tj<9r{x £ T : u(x) = a or u{x) = b} < +oo (4.1) 

(the number of points on the boundary in the topology of T of the active set is finite). 
Then the optimal control u belongs to H l / 2+t (T) and y £ H 1+t (Pl). 

Proof. The proof follows the lines of that of Theorem [TTJ Since ip reg £ H 2+t * (fl) and 
(p reg = 0 on T its normal derivative will be in H 1 ^ 2+t (T) provided t < t* and t < 1 (the 
condition (p reg = 0 on T is needed for t* > 1/2 to prove that the normal derivative tends 
to zero at the corners, and hence it is continuous; notice that for t* = 1 this continuity 
is not enough to have that the normal derivative is in i/ 3 / 2 (r).) 

This H l / 2+t (T) regularity is not affected by the projection formula (13.11) because t < 1 
and we are supposing m- The same happens with the singular terms such that cy i ^ 0 
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and Xj < 1. The rest of the singular terms in the expression for the normal derivative of 
the adjoint state will be in if 1 / 2+t (T) since t < m\j — 1 (see equations (13.81) and (13.91) 
for the expression of the normal derivatives of the singular part.) 

Let us prove that y £ H 1+ *(Q). Since u € H l / 2+t { T), there exists some U £ H 1+t (Q) 
such that U = u on T. Moreover, A U £ H~ 1+t (Q). So z = y — U is the solution of the 
boundary value problem 


—Az = A U in hi, z = 0 on T. 

Using the regularity results in [7J §23.C] (see also JIT], Th. 2.4.3 and §2.7]) we have that 
there exist a unique z reg £ H 1+t (Ll) and unique coefficients cy m such that 

3 

z — z reg T ^ ^ ^ .mLj sm(mXjdj). 

m=1 j£$™ 1+t 

Since t < X\, the singular part is in H 1+t (Ll), and so is the optimal state. □ 

Finally, we will describe the adjoint state for even more regular data. In the rest of 
this section we will suppose yn, £ IF 1,P *(U), p* > 2. 

For p > 1 and m £ Z define 

J ™ p = |j £ {1,..., M} such that 0 < mXj <3-and mXj 0 Z 

Now we have that J ™ = 0 if m > 5. We have to add the condition mXj £ Z otherwise 
logarithmic terms may appear. Define also: 

f 2 

L= < j £ {1,... , M} such that 0 < mXj <3-and mXj £ Z 

A direct calculation gives us that = 0 if m = 2 or m > 4, C {j : ujj = tt/2} 

and hf p c{j: ujj = 3ir/2}, and hence mXj = 2 if j £ L"' p . 

Consider now p > 2 such that, for m = 1,2,3 

p < p *, p < pd, P < 77 ——~r ^ 3 € V s* < oo. (4.2) 

2 — mXj 

In addition, we need to to assume 

^Z Vj £ {1,..., M}. 

XjP 

With this notation, we have the following result. 
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Corollary 17. There exist a unique function ip reg £ IF 3,p (n) and unique real numbers 
Cj jTn and dj,m such that 

5 

T = Preg + X X Omipf 3 sia(m\j6j) 
m=x 

+ E E dj^fj'Tj (log (rj) sin(2 6j) + 9j cos(2 9j)). 

m=l,3j£L™ p 


Proof. Since p* >2, yn £ L s * (Ll) for any s* < oo, and hence we can apply Theorem fill 
and we have that also y £ W 1,P (D). The result follows directly from the adjoint state 
equation (|3.3D thanks to the regularity result [10, Th. 5.1.3.5] 

5 

V = Treg + X X ^'X' sin ( mA J%) 
m=l jel^p 

+ X X dj^jr™^ (log(rj)sm(m\j9j) + 9j cos(m\j6j)) . 

m=l,3j£h™ v 

and using that m\j = 2 if j £ L™ . □ 

5 Problems without control constraints 

For problems without control constraints, we obtain similar results. Indeed, for convex 
domains and data yn £ L s (LI), s* > 2, it is obvious from Theorem fl8l below that the 
optimal control is a bounded function and hence all the results stated before apply. 
Nevertheless, for nonconvex domains, we will not obtain a continuous control and the 
singularities must be taken into account near all the corners. Therefore, the indexes for 
the expansion of the singular parts must be taken running through all the sets J™, and 
not only through the sets H™. 

Theorem 18. Suppose now that —a = b = oo and yn £ L s (Ll), s* > 2. Then there 
exists a unique (p reg £ W 2,s (Ll), u reg £ VF 1_1//s,s (r), y reg £ W 1,s (fi), for all s < s*, 
s <pd, and unique real numbers and (cj)a ; <i_ 2 /s such that 


3 


<P = Preg + X X 

(5.1) 

m=l jeJ™ 



u(x ) — U reg + 'y ] ''y ] aj t mf,jT 3 

+ X] Xj a j,m^j1" 3 

(5.2) 

m=l,3 jeJ™ 

tell 


3 

Dreg + X] X] a j,m^j r j 3 s j,m(9j) 

+ X OV/' sin ( A i^) 

(5.3) 

m =i jely 1 

\j<l—2/s 
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where ap m and Sj^ m (6) are given by the formulas (13.1311 and ()3.14l) . 

Proof. The proof is very similar to those of theorems [9l [FT] and [13j We will only empha¬ 
size on the main difference: at the beginning of the proof of Theorem [9] we used that the 
optimal control was bounded to obtain that the optimal state was a function in L s * (12). 
Now the optimal control is not bounded, so we use a bootstrapping argument to show 
that y G L S (I2) for s < s*, s < pn- The result follows then using the same techniques as 
before. Notice that now we do not need the sets H™, since we do not have to exclude in 
the expression of the singular part of the control the corners where the normal derivative 
of the adjoint state is not bounded. 

In a first step we have that u G L 2 (T), and hence y G H 1//2 (f2) C L 4 (I2). If s* < 4 the 
proof is complete since 4 < pn for any polygonal domain. 

Suppose s* > 4. The normal derivative of the regular part of the adjoint state is 
in W 1_1 / 4 ’ 4 (r), but now, since we have no control constraints, we have to take into 
account the normal derivative of the singular part near the non-convex corners. We 
have so far that the optimal control can be written as the sum of a regular part, which is 
in ty 1 - 1 /M(r) pi us a singular part, that behaves as r^ 1-1 . For the regular part we apply 
Lemma [3] and for the singular part we apply Lemma [5j and we have that the optimal state 
can be written as the sum of a regular part which is in IT 1,4 (II) C L s * (II) plus a singular 
part that behaves at worst as G L s {Pt) for all s < 2/(1 — min{l, Ai}) = po- So 

we have that y G L s (0) for all s < s*, s < pn- □ 

We will finish this section stating some regularity results of the optimal solution in 
the unconstrained case in some special situations. 

Corollary 19. Suppose the assumptions of Theorem [73 are satisfied. Then, for all 
p<s*,p<p n we have that (p G W 2 ’ P (D), u G W i-Vw^T) and y G W^ p {p). 

Proof. Since p < s* and p < pn < Pd, the regular parts of the involved functions satisfy 
u reg G VF 1-1 / p ’ p (T), y reg G W 1,P (Q) and (p reg G W 2,S (D) due to Theorem fl8l 

On the other hand, the assumption p < pn implies £ir^ 1-1 G W 1,p (Ll), and hence 
obviously £irjh G W 2,p (Ll) and £ir^ 1-1 G W l ~ l / p,p (T). Since these are the worst terms 
we may find in the singular parts, the proof is complete. □ 

With the same techniques of Section 21 (Corollary [TB] and [TUI) we can obtain the 
following result. 

Corollary 20. Suppose the assumptions of Theorem [ 73 are satisfied and yn G H 1 (12) 
for some t* < 1. Define t > 0 such that 

t <t*,t < 1, t < Ai — 1. 

Then u G H l / 2+t {T) and y G H l+t (i 2). 
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